Numerical Simulation of Impact Loading on Solid Propellant Within an Orthotropic Shell

The strain–stress state of the solid propellant rocket engines (SPREs) is simulated under impact. The effect of orientation of elastic and strength properties of orthotropic organoplastic shell material on the strain–stress state of the solid propellant is investigated. Normal and oblique impact of single steel cylinder projectiles, both simultaneous and at different times of multiple, converging steel spheric particles with SPRE are investigated in this study. The investigation is conducted numerically. The numerical modeling was carried out in a three-dimensional formulation by the method of finite elements for the continuous approach of the mechanics of a deformable solid. The destruction of the anisotropic material is described by the tensor-polynomial criterion of the fourth degree, which takes into account the influence of hydrostatic pressure.     

脉冲爆震火箭发动机引射模态的起爆实验研究

为实现脉冲爆震火箭发动机(PDRE)引射模态下主爆震室的起爆,采用航空煤油和氧气作为推进剂,设计了PDRE引射模态的模型机,采用压电传感器测量主爆震室中爆震波的压力和速度.在主爆震室中成功实现了5 ～8Hz稳定连续的爆震,爆震波的峰值压力能够达到3MPa,爆震波以1600～ 2000m/s左右的速度在主爆震室中传播.实验结果表明:PDRE引射模态下主爆震室的DDT距离,远低于常规高能电喷起爆下的两相PDRE的DDT距离;高频PDRE引射模态下主爆震室的起爆难度加大;加长主爆震室、末端增加收敛段可以提高引射模态的爆震性能.     

凹面腔内二维激波会聚特性研究

为了研究凹面腔内二维激波会聚的特性,设计、建立了二维激波会聚冷态实验系统,进行了二维激波会聚冷态实验,结果表明:当射流入口宽度增大时,相同射流入射压力下,喷入凹面腔内的射流强度增大,导致凹面腔内的激波会聚强度增强,从而压力脉动幅值增大,频率也随之增大.数值模拟了二维冷态激波会聚过程,从模拟结果中的压力和频率上看,模拟结果与实验结果基本上是吻合的.     

激波聚焦诱导气液两相爆震燃烧的数值模拟

对以激波聚焦和增加障碍物方式诱导煤油-空气气液两相爆震燃烧的过程进行了数值模拟.采用欧拉-拉格朗日方法建立了脉冲爆震发动机(PDE)中气液两相流的喷射、雾化、掺混过程.研究发现环形爆震波在爆震管凹腔内经过反射、汇聚后能够引燃可燃混合物.而在障碍物处,激波的反射和再反射聚焦能够形成高温高压点(2700K,25MPa),产生局部爆炸,有助于形成稳定的脉冲爆震燃烧(波面速度为1900m/s,温度为2 950K),有效地缩短由缓燃向爆震转变(DDT)距离至0.45m.     

喷嘴位置对脉冲爆轰发动机性能影响的试验研究

为了研究喷嘴位置对脉冲爆轰发动机性能的影响,设计加工了直径为80mm的汽油/空气两相脉冲爆轰发动机（PDE）,在PDE文氏管内不同位置安装喷嘴,进行了冷态雾化和热态燃烧转爆轰试验,分析了马尔文激光粒度仪测得的液滴尺寸分布和动态压力传感器测得的信号。试验结果表明：喷嘴的安装位置对脉冲爆轰发动机的性能具有显著影响;喷嘴安装在文氏管喉部时爆轰管内雾化混合效果最好,爆轰波峰值压力最大为3．5 MPa,工作频率最高为30Hz;喷嘴安装在文氏管入口时爆轰管内雾化混合效果最差,爆轰波峰值压力最小为0．9MPa,工作频率最低为15 Hz;在试验范围内,改善雾化混合效果有利于提高脉冲爆轰发动机的工作频率。研究结果对脉冲爆轰发动机的设计具有参考价值。     

气相爆轰波传播特性的数值模拟及实验对照

本文应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波进行了二维数值模拟.氢氧混合物的化学反应模型考虑了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.本文首先对爆轰波数值结果和实验结果进行了对照验证,然后对爆轰波在楔面反射由马赫反射向规则反射转变的过程进行了数值分析,得到了反射转变临界角,并和实验结果及理论分析进行了比较,结果是令人满意的;本文还对爆轰波的多波结构进行了初步的数值分析.     

脉冲爆震发动机模型机平均推力特性计算与实验研究

针对如何确定脉冲爆震发动机平均推力的问题,介绍了体积比冲估算法、推力壁压力曲线积分法和弹簧-质量-阻尼系统法的计算与测量原理。建立了PDE模型机推力测试台架系统。同时采用推力壁压力曲线积分法和弹簧-质量-阻尼系统法对一PDE模型机在不同频率多循环工况下的平均推力特性进行实验研究。对比分析了体积比冲估算法、推力壁压力曲线积分法和弹簧-质量-阻尼系统法所获得的PDE模型机平均推力计算值与实验值,发现：PDE模型机的平均推力随工作频率的升高基本上呈线性增加;同一工作频率下,采用体积比冲估算法得到的平均推力最大,推力壁压力曲线积分法次之,弹簧-质量-阻尼系统法最小;PDE模型机平均推力应处于推力壁压力曲线积分法与弹簧-质量-阻尼系统法所测推力之间。     

斜爆轰波的波角和法向速度-曲率关系初探

为研究斜劈诱导斜爆轰波的波阵面弯曲效应，以期为斜爆轰的不稳定性及其演化规律提供新的见解，基于加权本质无振荡（WENO）格式空间离散和附加Runge-Kutta方法时间离散的求解器，针对不同的化学反应参数（释热量、放热速率和化学反应区参考长度）条件，开展斜爆轰波的数值计算研究。结果表明斜爆轰波沿波阵面的波角变化可分为3个区域：区域I，波角平滑减小；区域II，波角跃升后衰减；区域III，波角有规律振荡。波阵面法向速度-曲率关系在区域I呈现准垂直直线变化趋势，并伴随着爆轰波强度的不断衰减；在区域III则呈现出"D"形曲线，即由极曲线段、光滑水平变化段和拟线性变化段组成，为类胞格结构的周期性演变；区域II可认为是以上两个区域特征的耦合。不同的化学反应参数对斜爆轰波波阵面的弯曲效应影响存在较大差别。     

燃料分布对旋转爆震波传播特性影响

为定量研究燃料分布对旋转爆震波传播特性的影响，针对两种环缝/小孔喷注方案开展了数值和实验研究。针对燃料(H2)和氧化剂(air)分开喷注的旋转爆震燃烧室模型，开展冷流掺混的三维数值模拟研究，给出了掺混均匀度评价参数，获得了两种喷注模型的反应物掺混均匀度沿轴向的变化规律。针对这两种喷注模型，开展旋转爆震波传播特性实验研究，分析了旋转爆震波的传播速度、工作稳定性、当量比边界等。研究结果表明:将燃料喷注小孔前移，可大幅提高燃烧室头部的掺混均匀度；随着掺混均匀度的提高，旋转爆震波的传播速度增加，传播稳定性明显提高，稳定工作的当量比下限从1.08拓宽至0.57。研究结论可为旋转爆震发动机喷注结构设计提供参考。      

曲面激波诱导斜爆轰的数值模拟

为了研究恰当比预混氢气-空气斜爆轰流场的波系结构和流动特征，基于带化学反应的Navier-Stokes方程，对弹头及楔穿越预混气体时诱导的斜爆轰进行了数值模拟。对流项的离散采用Steger-Warming格式，时间项采用二阶Runge-kutta方法。结果表明，对于弹头：（1）在亚爆轰条件下，能够模拟氢气-空气预混诱导爆轰流场的精细结构；（2）在超爆轰条件下，通过精细调整网格，能够很好地分辨强烈耦合的激波和燃烧波，且与Lehr实验吻合良好。对于楔结构：捕获到了清晰的三波点及其复杂精细的斜爆轰流场结构，预测的诱导燃烧距离、激波角和斜爆轰角与实验吻合良好；通过对流场波系结构变化过程的研究，获得了流场三波点随时间的演化过程。